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Abstract. Large f/-selfadjoint random matrices are considered. The matrix 
H is assumed to have one negative eigenvalue, hence the matrix in question has 
precisely one eigenvalue of nonpositive type. It is showed that this eigenvalue 
converges in probability to a deterministic limit. The weak limit of distribution 
of the real eigenvalues is investigated as well. 



Introduction 

The main object of this survey are random matrices that are not symmetric, 
but are selfadjoint with respect to an indefinite inner product. Spectrum and 
numerical range of some classes of such matrices were considered recently in [23ll24| . 
although, in the present paper we consider different instances. The inefinite linear 
algebra motivation of the present research is the follwoing. Consider an invertible, 
hermitian-symmetric matrix H e C nXTl . We say that X £ c™ x « j s 77-self adjoint 
if X* H = HX. This is the same as to say that A is selfadjoint with respect to an 
inner product 

[x,y]n ■= y*Hx, x,y e C™. 

Note that this inner product is not positive definite if H has negative eigenvalues. 
In the literature the space C" with the inner product [-,-]h is also called II K -space 
(where k is the number of negative eigenvalues of H) or Pontryagin space, the 
infinite dimensional case is considered as well [5J HO] . In the present paper the case 
when 

" -1 
In 

is considered. It is easy to check that for such H each i/-selfadjoint matrix has the 
form 

" a -b* 
b C 



(1) H = 



(2) X 



with x € K, b E C N and a hermitian-symmetric matrix C G C NxN . Due to the 
famous theorem of Pontryagin |27j the matrix X has precisely one eigenvalue /3, 
for which the corresponding eigenvector x satisfies [a;,x]jy < 0. The problem of 
tracking the nonpositive eigenvalue was considered for example in [9, 29 . In those 
papers the setting was non-random and X was in the family of one dimensional 
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Figure 1 . Eigenvalues of a random matrix X iqo computed with R [21] 

extensions of a fixed operator in an infinite dimensional Ili-space. The aim of the 
present work is to investigate the behavior of /3 when X is a large random matrix. 
We show that the main method of [9J [29] - the use of Nevanlinna functions with 
one negative square - can be adapted to the random setting as well. 

A classical result of Wigner [3T] says that if the random variables y^, < i < 
j < +oo are real, i.i.d with mean zero and variance equal one, then the distribution 
of eigenvalues of a matrix 

Y N = -f=[Vii]y=o> 

where yji = for j > i, converge weakly in probability to the Wigner semi- 
circle measure. Note that by multiplying the first row of Yjv by -1 we obtain a 
-ff-selfadjoint matrix Xn- A result of a preliminary numerical experiment with 
gaussian y^ is plotted in Figure [T] Note that the spectrum of Xn is real, except 
two eigenvalues, lying symmetrically with respect to the real line. Although we pay 
a special attention to the above case, we study the behavior of the eigenvalue of 
nonpositive type in a more general setting. Namely, we assume that the random 
matrix Xn in <C NxN is of a form 

„ _ a/v — b* N 
I bjy Cn 'J 

with oat, bpf and Cn being independent. Furthermore, the vector 6jy is a column 
of a Wigner matrix and ajy converges weakly to zero. The only assumption on Cn 
is that the limit distribution of its eigenvalues converge weakly in probability, see 
(R0)-(R3) for details. In Theorem [6] we prove that under these assumptions the 
non-real eigenvlaues converge in probability to a deterministic limit that can be 
computed knowing the limit distribution of eigenvalues of Cn- In the case when 
Cn is a Wiegner matrix the nonreal eigenvalues converge to ±iV2/2, cf. Theorem 
[8] Furthe rmore, under a technical assumption of continuity of the entries of Xn, 
we show in Theorem [7] that the limit distribution of the real eigenvalues of Xn 
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coincides with the limit distribution of eigenvalues for the matrices Cjv- Again, in 
the case when Cjv is a Wiegner matrix we obtain a more precise result. Namely, in 
Theorems we show that the real eigenvalues . . . , of Xn and the eigenvalues 
, . . . , Ajy of C'n satisfy the following inequalities: 

It shows that the nonreal eigenvalue of Xn plays an analogue role as the largest 
eigenvalue in one-dimensional, symmetric perturbations of Wiegner matrices. This 
fact relates the present paper to the current work on finite dimensional perturba- 
tions of random matrices, see [21 El SI El EJ [T3J [TS] and references therein. Also note 
that Xjv is a product of a random and deterministic matrix, such products were 
already considered in the literature, see e.g. [3D]. 

The author is indebted to Anna Szczepanek for preliminary numerical simula- 
tions and Maxim Derevyagin for his valuable comments. Special thanks to Anna 
Kula, Patryk Pagacz and Janusz Wysoczahski. 

1. Functions of class A/i 

The Nevanlinna functions with negative squares play a similar role for the class 
of iZ-self adjoint matrices as the class of ordinary Nevanlinna function plays for 
hermitian-symmetric matrices. This phenomenon has its roots in in operator the- 
ory, we refer the reader to [5] HI HFJ El US] and papers quoted therein for a precise 
description of a relation between A/" K -functions and selfadjoint operators in Krein 
and Pontryagin spaces. We begin with a very general definition of the class A/"«, 
but we immediately restrict ourselves to certain subclasses of those functions. 

We say that Q is a generalized Nevanlinna function of class Af K \16\ 1 1 9 j if it is 
meromorphic in the upper half-plane C + and the kernel 



Q(z) - Q{w) 

N(z,w) = z 

z — w 

has precisely n negative squares, that is for any finite sequence z\, . . . , Zk € C + the 
hermitian-symmetric matrix 

has not more then k nonpositive eigenvalues and for some choice of zi, . . . , z^ it 
has precisely k nonpositive eigenvalues. In the present paper we use this definition 
with k = 0,1. 

The class A/o is the class of ordinary Nevanlinna functions, i.e. the functions that 
are holomorphic in C + with nonnegative imaginary part. By M^(M) we denote the 
set of positive, bounded Borel measures on K. For fx <E M^" (R) we define the 
Stieltjes transform as 



fi(z) = j dt, zeC\supp/i. 



Clearly, jj, belongs to the class A/o and the values of fj, in the upper half-plane 
determine the measure uniquelly by the Stieltjes inversion formula. Although not 
every function of class A/q is a Stieltjes transform of a Borel measure (cf. [H]), this 
subclass of A/o functions will be sufficient for present reasonings. Also, we will be 
interested in a special subclass of A/i functions, namely in the functions of the form 
([3]) below. We refer the reader to the literature [5] [TDJ for representation theorems 
for Af K functions. 
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Proposition 1. If fi € A/ b + (M), a e K then 

(3) Q(z) = fh(z) + a — z 

is a holomorphic function in C + and belongs to the class M\. Furthermore, there 
exists precisely one z € C such that either z € C + and 

(4) Q(«,) = 0, 
or zo € K and 

(5) lim e (-oo,0]. 

2^-20 2 — ZQ 

The symbol A above denotes the non-tangential limit: 

zeC+, z z , tt/2 - 6< < arg(z-z ) < vr/2 + 0, 

with some € (0,7r/2). We call zo € C + U K the generalized zero of nonpositive 
type (GZNT) of Q(z). The first part of the Proposition can be found e.g. in [17] . 
while for the proof of the 'Furthermore' part in the general context)^] we refer the 
reader to [191 Theorem 3.1, Theorem 3.1']. In view of the above proposition we can 
define a function 

G : M+(R) x!^C + 

by saying that G(/i, a) is the GZNT of the function fi(z) + a — z. The following 
proposition plays a crucial role in our arguments. 

Proposition 2. The function G is jointly continuous with respect to the weak 
topology on M^(IR) and the standard topology on K. 

Proof. Assume that (/i n )„ C (M) converges weakly to /i g M^(R) and n„el 
converges to a £ R with n — > oo. Take a compact if in the open upper half-plane, 
with nonempty interior. Then fj, n converges uniformly to fi on the set K . Indeed, 
if r = sup teR zeK l/\t — z\ then 

sup \jj, n (z) - Ao0)| < r\fi n - /i |(R), 
zeK 

the latter clearly converging to zero with n —¥ oo. In consequence, £i n (z) + a n — z 
converges to fi(z) + a — z uniformly on K with n — > oo. By [3T] the GZNT of 
An( z ) + a n — z converges to the GZNT of fl(z) + a — z, which finishes the proof. □ 

2. ii— SELFADJOINT MATRICES 

In this section we review basic properties of selfadjoint matrices in indefinite 
inner product spaces introducing the concept of a canonical form and showing its 
relation with Wi -functions. Let H G c( n+1 > x (™ +1 ) (n € N \ {0}) be an invertible, 
Hermitian-symmetric matrix. We say that X 6 c(n+i)x(n+i) j s jj -selfadjoint if 
X*H = HX. Our main interest will lie in the matrix 



(6) H 



-1 

o /„ 



^For arbitrary A/i function zo = oo can be also the GZNT, in that case Hindoo zQ(z) S [0, oo). 
However, this is clearly not possible for Q of the form ([31. 
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where /„ denotes the identity matrix of size n x n. As it was already mentioned, 
each H-self adjoint matrix has the form 

a -b* 
b C 



(7) 



X 



with a £ M, b £ C™ and hermitian-symmetric C £ C nxn . Due to [T3] there exists 
an invertible matrix S and a pair of matrices H',S' £ c("+ 1 ) x (™+ 1 ) such that 
X = S^X'S H = S*H'S and A', H' are of one of the following forms: 
Case 1. 



Case 2. 



Case 3. 



Case 4. 



X' = 
with /3 £C 

with ft £ 

X' = 
with ft £ 



ft 

p 
, Ca,. 



X' -- 

a, Ci 



©diag(C 2 , 
■ , Cn e a. 

®diag(Ci, 

, C™ e a. 



■ , Cn), # — 



,C„), ff' = [-!]©/„ 



p I 

/3 

C2, ■ ■ 



Bdiag(C 2 , ..-,Cn), 

,C«eR,7e {-1,1} 



H' 



X' 



diag(C 3) ...,Cn), H 1 



In — ? 



With/3GK, C3,---,Cn € 



The pair is 
14"] for the proof and 



It is easy to verify that in each case X' is -ff'-symmetric. 
called the canonical form of (X,H). We refer the reader to 
for canonical forms for general H-symmetric matrices and to [50] for the infinite- 
dimensional counterpart of the theory. At this point is enough to mention that the 
canonical form is uniquely determined (up to permutations of the numbers Q) for 
each pair (X, H), where X is if-selfadjoint. Note that in each of the cases P is an 
eigenvalue of X and there exists a corresponding eigenvector x £ C n+1 satisfying 
[a;, x]h < 0, furthermore, P is the only eigenvalue in C + U R having this property. 
Therefore, we will call ft the eigenvalue of nonpositive type of X. 
Observe that the function 

(8) Q(z) = a-z + b*(C-z)- 1 b 

is an A/i -function. Indeed, if D = UCU* = diag(Ai, . . . , A„) is a diagonalization of 
the hermitian-symmetric matrix C and d = Ub then 



Q(z) = a — z + — — — = a — z + fi(z), where /1 = \dj 



A, 

=1 3 



3=1 " 3 = 1 

and we may apply Proposition[l] The following lemma is a standard in the indefinite 
linear algebra theory. We present the proof for the reader's convenience. 

Lemma 3. Let X andQ be defined by (JtJ) and respectively. A point ft £ C + UM 
is the eigenvalue of nonpositive type of X if and only if it is the GZNT of Q(z). 
Furthermore, the algebraic multiplicity of ft as an eigenvalue of X equals the order 
of p as a zero ofQ(z). 
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Proof. First note that, due to the Shur complement formul^ 

where e denotes the first vector of the canonical basis of C™ +1 . Let (X',H') be 
the canonical form of (X, H) and let S be the appropriate transformation. Conse- 
quently, 

(9) --^ ) ={Se l )*H'{X'-z)- l Se 1 . 

Below we evaluate this expression in each of the Cases 1-4. Let / = [/o, . . . , f n ] T — 
Se. Note that 

(10) KH'h = elHex = -1, 
independently on the Case. 



Case 1. Observe that /o/i 7^ 0, otherwise f*H f > 0, which contradicts (10 1. 
Due to (|9| one has 

1 _ /0/1 , /1/0 , y-^ \fj\ 2 



Q(z) p-z p-z f^Q-z 

Hence, f3 £ C + is a simple pole of — 1/Q and consequently it is the GZNT of Q and 
a simple zero of Q. 

Case 2. Observe that |/o| 2 > 2j=i \ fj\ 2 > otherwise f*H'f > 0, which contradicts 

( [Tot ). Due t0 ® one nas 

1 -l/ol 2 , y LM! 

Hence, the residue of — 1/Q in /3 is less then zero. Consequently Q(f3) = 0, Q'(/3) < 
and is the GZNT of Q. 

Case 3. Observe that |/i| 2 > 0, otherwise f*H'f > 0, which contradicts (10 1. 
Due to Q one has 



2 7 Re/o/! , -7I/1I 2 , J2 l^l 2 



z 



Hence, P is pole of -1/Q of order 2. Consequently, Q(/3) = Q'{fi) = 0, Q"(^) ^ 
and /3 is the GZNT. 

Case 4. Observe that |/ 2 | 2 > 0, otherwise f*H'f > 0, which contradicts (10). 
Due to ^ one has 

1 _2Re/ Q / 2 + |/ 1 | 2 , -2ReA/ 2 , |/ 2 | 2 " |/,| 2 



Cfl _ ?\2 + (R_ ?\3 + 2-^ 



Q(z) p-z {p-z)* {fi-zf f^Q-z' 

Hence, P is pole of -1/Q of order 3. Consequently, Q(P) = = Q"(/3) = 0, 

Q'"C9) + and /? is the GZNT of Q. □ 



2 It is well known |17| that —1/Q belongs to A/i provided that Q belongs to JVi , however, this 
information is not essential for the proof. 
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3. Random .H-selfadjoint matrices 



By Xjf, Hn we understand the following pair of a random and deterministic 
matrix in c(n+i)*(n+i) 



(11) 



Xjsi — 



a N 
b N 



N 



H]y = 



-1 







In 



where a at is a real- valued random variable, 
a hermitian-symmetric random matrix in d 



6 at is a random vector in 
: nxn_ Note that x 



< \ N 



N , and Cat is 
N is i?AT-symmctric. 



By Af < 
random measure on 

Recall that 
(12) 



we denote the eigenvalues of Cat and by Vn we denote the 



v N 



N 



v N {z) 



— > A' 



tr(CW - z)- 
7V 



The assumptions on Xjy are as follows: 

(RO) The random variable a at is independent on the entries of the vector and 

on the entries of the matrix Cat for each N > 0, furthermore oln converges 

with N — > oo to zero in probability. 
(Rl) The random vector 6^ is of the form 

bN '■= —7=[ x io]i=t,...,N, 



N 

where [a;jo]j>o are i-i-d. random variables, independent on the entries of 
Cn for N > 0, of zero mean with £7|cCjo | 2 = s 2 for j > 0. 
(R2) The random measure vn converges with N — > oo to some non-random 
measure /io weakly in probability 

All the results below hold also in the case when all variables Xjo (j > 0) are real, 
in this situation b* N is just the transpose of 6at- The entries of Cn might be as well 
real or complex. In Section [5] we will consider two instances of the matrix Cn- a 
Wiegner matrix and a diagonal matrix. In the case when Cn is a Wiegner matrix 
the proposition below is a consequence of the isotropic semicircle law [T2 US] . We 
present below a simple proof of the general case, based on the ideas in |22j . 

Proposition 4. Assume that (Rl) and (R2) are satisfied. Then for each z G C + 



N 



l b 



N 



flo(z) (iV-»oo) 



in probability. 



By 1 1 y 1 1 we denote the euclidean norm of y G 



Proof. In the light of Chebyshev's inequality, ( 12 ) and assumption (R2) it is enough 
to show that 

2 



(13) 

Observe that 
(14) 



lim E 

N~yoo 



b* N {C N -z)- x b N 



2 tr{C N -z)- 1 



N 



0. 



E 



b* N (C N - z) l b N - s : 



,ti(C N - z)- 1 
N 
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E\b* N (C N - z) _1 &jv] - S 4 £ 



tr(C N -z)- 1 



N 



(15) 



2 Re E Is 



2 tr(C N -z) 1 ( h * lr 2 tr(CAr-z) 1 

o n {Cn — z) On — s - 



N 



N 



First we prove that the summand (15) equals zero. Indeed, conditioning on the 



cr-algebra generated by the entries of the matrix Cjv and setting 



[cij\? j=1 = (C N - z)- 1 



one obtains 



N 



N 



JV 



E KEf E 



iv 



E 



iV 



JV / JV 2 



N 



^KE^(E^4- s2 E 



N , 

»=i \j=i 



3=1 



iV 



Next, observe that 



N 



E\b* N (C N - zy^N] 2 =e Y, 

ijkl—1 



x 0i x 0j x 0k x 0l 
N 2 



4 \ E(cgCjj) 4 \ - E(cjjCjj) _ 4 
S 2^ 7V2 + S TV 2 _ 

This allows us to estimate (fl4|) by 

iV 



tr(C. 



JV 



N 



/v 



^E 



C ij C ij 

N 2 



*E 

ij=i 



Cj t Ci 



N 2 



< s E- 



(C N - z)- 



N 



s 4 dist(z, er(CV)) 2 s 4 



N 



Im z 2 N ' 



which finishes the proof of (13). 



□ 



Let Um be a unitary matrix, such that UnCnU n is diagonal and let aV 
[d^, . . . , dj^] T = f/jv&jv- Denote by {In the measure defined by 

JV 

mjv = 

3=1 



and observe that /ijv(^) = b* N (CN — z) b 



Proposition 5. Assume that (Rl) and (R2) are satisfied. Then the sequence of 
random measures [an converge weakly with N — > oo to fiQ in probability. 
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Proof. First note that almost surely /zjv(R) — > s fi (W.) with N — > oo. Indeed, 

iV AT 

mm) - E Ki 2 = ii^ii 2 = im 2 = n E i^-i 2 ' 

J'=l .7=1 

which converges almost surely to s 2 by the strong law of large numbers. Further- 
more, Proposition [4] shows that /2jv(z) converges in probability to flo(z) for every 
z G C + . Repeating the proof of Theorem 2.4.4 of [T] we get the weak convegence 
of \x 7v in probability. □ 

4. Main results 

Theorem 6. If (RO) - (R2) are satisfied then the eigenvalue of nonpositive type 
/3jv of Xn converges in probability to the GZNT /3 of the N\ -function 

Qo{z) = -z + s 2 fi (z). 

Proof. Consider a sequence of Afi -functions 

(16) Qn(z) = fljv - z + fi N (z). 

Recall that each of those functions has precisely one GZNT which, by definition 
of /ijv and Lemma |3j is the eigenvalue of nonpositive type /3/v of X^. Recall that 
ajv converges to zero in probability by (RO) and /ijv converges to fj,Q in probabil- 
ity by Proposition [3J Let d be any metric that metrizises the topology of weak 
convergence on M^(M). Since (3n is a continuous function of /z/v and ajy (Propo- 
sition [2]), for each e > one can find S > such that for each N > the event 
{\on\ < 6, d(fj,N,fj,o) < S} is contained in {\(3n — Po\ < £}• Using the assumed in 
(RO) independence of [In and one obtains 

P(\(3o ~ Pn\ >£)< P(\on\ > 5) ■ PidifiN^o) > S). 
Hence, /3jv converges to /3q in probability. □ 

As it was explained in Section]!] each matrix Xn has, besides the eigenvalue /3/v 
of nonpositive type, a set of real eigenvaules Q^ N , . . . , £jy , where fcjv = 1 in Case 1 
and 3, fc/v = 2 in Case 2 and kpj = 3 in Case 3. By Tpj we denote the empirical 
measure connected with these eigenvalues: 

^E€- 

j=k N 

Theorem 7. If (R0)-(R2) are saisfied and the random variables {xQj : j > 0} are 
continuous, then the measure rjv converges weakly in probability to /io- 

Proof. We use the notations U^,d N and fi N from the previous section, let also Qn 
be given by (16). Note that the set {y € C N : (L/jvy)j — 0} is of Lebesgue measure 



zero. Hence, with probability one d^ 7^ for j = 1,...,N, N > 0. Therefore, 
/t/v(z) is a rational function almost surely with poles of order one in , . . . , Ajy. 
Furthermore, 

_ , , («»-*) nf.,(\"-^)+Ef.,Ki 2 n ) v,(^-=) 
Qn{ > = itf-W-) ■ 

In consequence, Qn bas exactly N + 1 zeros counting multiplicities, all of them 
different from A^, . . . , Ajy. Due to the Schur complement argument, each of those 
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zeros is an eigenvalue of the matrix Xn € (t^(n+i)^(n+i) ^ Furthermore, due to 
Lemma [3] the algebraic multiplicity of /3n as eigenvalue of Xn equals the order of 
/3jv as a zero of Qn- In consequence, the spectrum of Xn coincides with the zeros 
of Qn and Pn is the only zero of order possibly greater then one[^] 

On the other hand, the function /tjv is increasing on the real line with simple poles 
in Xi , . . . , Ajy. Hence, in each of the intervals (Af, (j = 1, ■ ■ ■ , N — 1) there 

is an odd number of zeros of Qn, counting multiplicities. Consequently, in each of 
the intervals (A^, Xj +1 ) (j = 1, . . . , N — 1) there is precisely one zero of Qn, except 
possibly one interval that contains three zeros of Qn- Out of these three zeros of 
Qn either one or two of them belong to the set {Ck N t ■ ■ ■ > Civ}; accordingly to the 
canonical form of Xn- Hence, in each of the intervals (Af, A^_i) (j = 1, ■ ■ ■ ,N —1) 
there is precisely one of the eigenvalues ^ except possibly one interval 

that contains two of the eigenvalues (J? N , . . . ,£/v- Consequently, the weak limit of 
tn in probability equals the weak limit of i>n ■ 

□ 

5. TWO INSTANCES 

In the present section we consider two instances of Cn '■ the Wigner matrix and 
the diagonal matrix. These both cases appear naturally as applications of main 
results. We refer the reader to [55] for a scheme joining both examples. 
Consider an 77-self adjoint real Wigner matrix 

1 

with Xij real, Xij = Xji (0 < i < j < oo), i.i.d., of zero mean and variance 



(17) := #A^i]£ =0 , 



equal to s , and let Hn be defined as in (111. Clearly Xn is -Hjv-selfadjoint and 
satisfies (R0)-(R2) with /i equal to the Wiegner semicircle measure a. The Stieltjes 
transform of the a equals 



— z + \/ z 2 — 4s 2 

< Z ) = ^2 



2s 2 

It is easy to check that /3q — ^si is a zero of Qo{z) = — z + s 2 a(z). Hence, flo is 
the GZNT of Qq and we have proved the first part of the theorem below. 



Theorem 8. Let Xn be defined by (17). Then 

(i) Pn converges in probability to /?o = ^ s ' l \ 

(ii) if, additionally, the random variables bij (0 < i < j < oo) are continuous, 
then the probability of an event that there are precisely N— 1 real eigenvalues 
(,2 < ■ ■ ■ < ^ of Xn and the inequalities 

(18) Af <C 2 7V <A^<--.<A^ 1 <C A r<A^. 
are satisfied, converges with N to 1. 

Proof, (ii) Assume that 

(19) |Av-A)| < ~s. 

Then the canonical form of (Xn,Hn) is as in Case 1. In consequence, there are 
exactly N — 1 real eigenvalues > ■ ■ ■ > Cn °f -^-N- Let us recall now the arguments 



^In other words: e is almost surely a cyclic vector of X? 
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FIGURE 2. The real and imaginary part of /3jv, with real, gaussian 
entries of and s 2 = 1, computed with R [2"5] , 



from proof of Theorem [7] The function /tjv is increasing on the real line with 
simple poles in A^, . . . , Ajy. In each of the intervals (Aj^, Xf + \) (j = 1, • • • , N — 1) 
there at least one of the eigenvalues Cis, • • • , Cjv • Consequently, each of the intervals 
(A^, Aj^) (j = 1, . . . , N — 1) contains precisely one of the eigenvalues , . . . , Cat • 
To finish the proof it is enough to note that by point (i) for every e > there exists 



Nq > such that for N > No the probability of ( 19 ) is greater then 1 — e. □ 



The numerical simulations of values of Re and Im /3/v can be found in Figure 
[2] Note that /3 lies in open upper half-plane and Q is satisfied. We provide 
now an example when /3q £ R and show that each number in [0, oo) can be the 
limit in ([5|. Let — 0, Xio (i = 1, 2, ,...,) be independent real variables of zero 
mean and variance s 2 and let Cjv = diag(ci, . . . ,cjv)j where the random variables 
{cj : j = 1, . . . } are i.i.d. and independent on (i = 1,2,,... ,). Furthermore, let 
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Figure 3. The imaginary part of /3jv- 



the law of Cj (which is simultaneously the limit measure /Xq) be given by a density 




: t e [-1,1] 

: t 6R\ [-1,1] ' 



An easy calculation shows that 

lim Mfl = 3 . 

Hence, 

lini -z+^AoM = _ 1 + 3s2 

z^O z 

and the function 

Qo(z) = -z + fio(z) = -z+ f pQ-dt 

has a GZNT at z = if s 2 < 1/3. Note that /3o = lies in the support of n$. 
The case s 2 = 1/3 is plotted in Figure [3] Only the imaginary part is displayed, 
since the numerical computation of the real part of /3jv might be not reliable in case 
/3jv G K- One may observe that the convergence of /3jv is worse in Figure [2] Also, 
the canonical form of (X N , H N ) changes with N, contrary to the case when H N X N 
is a Wigner matrix. In the case s 2 < 1/3 in numerical simulations point /3/v is real 
for all N. 
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